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We apply the Bogoliubov compensation principle to QCD. The non-trivial solution ol compen- 
sation equations for a spontaneous generation of the anomalous three-gluon interaction leads to 
the determination of parameters of the theory, including behavior of the gauge coupling a B (Q 2 ) 
without the Landau singularity, the gluon condensate Vi — 0.01 GeV 4 , mass of the lightest glueball 
Mo ~ 1500 MeV in satisfactory agreement with the phenomenological knowledge. The results 
strongly support the applicability of N.N. Bogoliubov compensation approach to gauge theories of 
the Standard Model. 
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I. INTRODUCTION 

We are now sure, that QCD is the genuine theory of 
strong interactions. In the perturbative region at high 
momenta the theory excellently describes the totality of 
data. However at low momenta the perturbative theory 
fails. Firstly it is seen from the well known momentum 
dependence of the running coupling. Let us show here 
the three loop expression for a s (fi) 
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where A is the QCD scale parameter and 
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Thus from (JTJ) we see, that for /i = A3 we have the pole 
(and the cut at the same point). 

At first such pole was disclosed in QED [lj, |2( and 
thus was called the Landau pole. The existence of the 
pole makes a theory internally contradictory. As for 
QED, L.D. Landau himself in the issue dedicated to Niels 
Bohr had first stated, that for a realistic number of 
the charged elementary fields the pole was situated far 
beyond the Planck mass and so it presumably could be re- 
moved by quantum gravitation effects. However in QCD 
pole is situated in the observable region of few hundreds 
MeV. As far as we know, there is no way to get rid of 
such pole in the framework of the perturbation theory. 
It is a general belief, that non-perturbative contributions 
somehow exclude the pole. For reviews of different pos- 
sibilities see e.g. [1,0]. 

In the present work we would demonstrate just how 
the pole in ([1} could be eliminated in the approach to 
non-perturbative effects in gauge theories, which was in- 
duced by_the famous N.N. Bogoliubov compensation ap- 
proach 
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For low momenta region we take expression (JT]) with num- 
ber of flavors Nf = 3 and take for normalization its value 
at mass of r-lepton. We have 



a s (M T = Yin MeV) 
From here we obtain 



0.32 ±0.05. 



A 3 = (345 ± 19) MeV . 
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II. THE COMPENSATION EQUATION 

For the beginning we consider pure gluon QCD with- 
out quarks. We start with Lagrangian with gauge group 
SU(3). That is we define the gauge sector to be color 
octet of gluons F". 



L = -- F a F a ■ 

4 A 1 * V" 
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where we use the standard notations. Let us consider 
a possibility of spontaneous generation of the following 
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effective interaction 

G b 

3~I ' ^ abc ^ t " y ^ up ' ^ 

which is usually called the anomalous three-gluon inter- 
action. 

Here notation ^ • f abc F pv F^ p F£ means correspond- 
ing non-local vertex in the momentum space 

(27r) 4 G fabc{g^(q P pk - p P qk) + g vp (k^pq - q^pk) + 
+9 P p,(p^qk - k v pq) + q^k v pp - k^p v q p ) x 
xF(p,q,k)5(p + q + k) + ... ; (7) 

where F(p, q, k) is a form-factor and p, /i, a; q, v, b; k, p, c 
are respectively incoming momenta, Lorentz indices and 
color indices of gluons. 

In accordance to the Bogoliubov approach [5j, |6[ in ap- 
plication to QFT Q we look for a non-trivial solution of a 
compensation equation, which is formulated on the basis 
of the Bogoliubov procedure add — subtract. Namely 
let us write down the initial expression ([5]) in the follow- 
ing form 

L = Lq + Li n t j 

1 C 
T. n = p a p a J ■ f h F a F h F c ■ (%) 

^0 — ^ r fj.u r pv ' g| J abc r p V r vp r pp j V°y 

G b 
Lint ~ g| f abc ^fiv ^vp ^PP ' 00 

Here notation - § • f abc F" v F* p F c pu is already ex- 
plained ([7]). 

We mean also that there are present four-gluon, five- 
gluon and six-gluon vertices according to expression for 
Lp V ©■ Note, that inclusion of total gluon term F^ v F pv 
in the new free Lagrangian ([8]) is performed in view of 
maintaining the gauge invariance of the approach. 

Effective interaction ^ is called anomalous three- 
gluon interaction. Our interaction constant G is to be 
defined by the subsequent studies. 

Let us consider expression ([S]) as the new free La- 
grangian Lq, whereas expression © as the new interac- 
tion Lagrangian Li nt . It is important to note, that we 
put into the new free Lagrangian the full quadratic in F 
term including gluon self-interaction, because we prefer 
to maintain gauge invariance of the approximation being 
used. Indeed, we shall use both four-gluon term from the 
last term in ((8]) and triple one from the last but one term 
of ([5]). Then compensation conditions (see for details Q) 
will consist in demand of full connected three-boson ver- 
tices of the structure ([7]), following from Lagrangian Lq, 
to be zero. This demand gives a non-linear equation for 
form-factor F. 

Such equations according to terminology of works 
0, @] are called compensation equations. In a study 
of these equations it is always evident the existence of a 



perturbative trivial solution (in our case G = 0), but, in 
general, a non-perturbative non-trivial solution may also 
exist. Just the quest of a non-trivial solution inspires 
the main interest in such problems. One can not succeed 
in finding an exact non-trivial solution in a realistic the- 
ory, therefore it is of great importance to choose an ad- 
equate approach, the first non-perturbative approxima- 
tion of which describes the main features of the problem. 
Improvement of a precision of results is to be achieved 
by corrections to the initial first approximation. 

Thus our task is to formulate the first approximation. 
Here the experience acquired in the course of performing 
works 0-[| could be helpful. Now in view of obtaining 
the first approximation we would make the following as- 
sumptions. 

1) In compensation equation we restrict ourselves by 
terms with loop numbers 0, 1. 

2) We reduce thus obtained non-linear compensation 
equation to a linear integral equation. It means that in 
loop terms only one vertex contains the form-factor, be- 
ing defined above, while other vertices are considered to 
be point-like. In diagram form equation for form-factor 
F is presented in Fig.l. Here four-leg vertex correspond 
to interaction of four bosons due to our effective three- 
gluon interaction. In our approximation we take here the 
vertex with interaction constant proportional to g G. 

3) We integrate by angular variables of the 4-dimensional 
Euclidean space. The necessary rules are presented in 
paper 0. 

Let us note that such approximation was previously 
used in works [8l-[l0| in the study of spontaneous genera- 
tion of effective Nambu - Jona-Lasinio interaction. It was 
shown in the works that the results agree with data with 
average accuracy ~ 10 — 15%. Thus we could hope for 
such accuracy in the present problem. Let us formulate 
compensation equations in this approximation. 

For free Lagrangian Lq full connected three-boson ver- 
tices with Lorentz structure ([Jj are to vanish. One can 
succeed in obtaining analytic solutions for the following 
set of momentum variables (see Fig.l): left-hand legs 
have momenta p and —p, and a right-hand leg has zero 
momenta. However in our approximation we need form- 
factor F also for non-zero values of this momentum. We 
look for a solution with the following simple dependence 
on all three variables 

F( P uP,P 3 ) = F(d±A±A). (10 ) 

Really, expression (|10p is symmetric and it turns to F(x) 
for P3 = 0, p\ = p\ = x. We consider the representa- 
tion (|10p to be the first approximation and we it would be 
advisable to take into account the corresponding correc- 
tion in forthcoming studies. We shall also discuss below 
some possible corrections due to this problem. 
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FIG. 1: Diagrams, describing the compensation equation. Lines correspond to gluons, black circles correspond to vertex ([7]), 
open circles to the same vertex with unity form-factor, open circles with four legs correspond to vertex simple point 

corresponds to usual perturbative vertex. 



At first let us present the expression for four-boson vertex 
V(p,m,\; q,n,a; fc,r,r; Z,s,tt) / , 

t ( 27r )4 =9G[f f (U{k,l;a,T,TT,X) - U(k,l; X,T,n,a) - U{l,k;a,ir,T, A) + 

U{l,k;X,TT,T,a) + U(p,q;n, X,<t,t) - U(p,q;r, A, cr, tt) - U(q,p;ir,cr, A,r) + U(q,p; r, ct, A,7r)) + 

/ a ™ J""* (t% i; ff, A, tt, r) - tf(|,p; a, tt, A, r) - tf(p, /; r, A, tt, a) + (7(Z,p; r, tt, A, a) + U(k, q; tt, r, a, A) - (11) 

U(q,k; tt, a, r, A) - U(k, q; A, r, a, tt) + [% fc; A, a, r, tt)) - f asn f amr (U(k,p; a,T, A, tt) - U{p,k; a, A, r, tt) + 

U(p,k;Tr,X,T,a) -U(k,p;ir,T,X,a) - U(l,q;T,n,a, A) + U(l, q; A, tt, cr, r) - U(q,l; X,a,ir, t) +U(q,l; T ,a,Tr,X))) ; 

C/(fc, Z;ct,t, 7T, r) = (k^lrg^x - k a hg^r + K hg aT - {kl)g aT g^x) x ^(fc, ^, -(fc + 0) ■ 

Here triad p, m, A etc means correspondingly incoming momentum, color index, Lorentz index of a gluon and F is 
the same form-factor as in expression ([7]). 

Now according to the rules being stated above we obtain the following equation for form-factor F(x), which 
corresponds to Fig.l. 



C 2 N I r Y 
F {x) =-—A F(y)ydy 



/ F(y)y 3 dy + -L / y 2 dj, + | 

2x z J Q ox J 6 



a; 
12 



2 '"^u + ^r Ws 



2/ 



GgW y y (5x-6y) 



o 



24; 



tt 2 \ix (a: - 2y) 



F(y)dy 



F(y)dy 
(3x-4y) 2 (2y-3x) 



(12) 



3» x 2 (x — 2y) ^ 



GgN { [ Y 3(x 2 - 2y 2 ) , 3(4?/ - 3x) 2 (a; 2 - 4^ + 2y 2 ) 



V^x 8(2y - a 

3a; 2 — 4xy — 6y z 
16y 2 



8x 2 (2y- a;) 2 



5y 2 — 12a;y 
16a; 2 



F(y)dy + 



F{y)dy 



V 



Here x = p 2 and y = q 2 , where q is an integration the initial equation by coupling constant G in view of 
momentum, number of colors N = 3. Here we also divide looking for non-trivial solutions. Of course, the trivial 
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solution G — is always possible. 

The last four terms in brackets represent diagrams with 
one usual gauge vertex (see three last diagrams at Fig.l) 
These terms maintain the gauge invariance of results 
in this approximation. Note that one can additionally 
check the gauge invariance by introduction of longitudi- 
nal term di k p k v / (k 2 ) 2 in boson propagators to verify in- 
dependence of results on di in this approximation. Ghost 
contributions also give zero result in the present approx- 
imation due to vertex ([7]) being transverse: 

P^V(p, 1, & W = q»V(p, q, k) Mvp = k p V(p, q, k)^ = 0; 
V(p, q, k) pvp = g pv (q p pk - p p qk) + g„ p (k p pq - (13) 
qnpk) + g PP {puqk - k v pq) + q p k v p p - k p p v q p . 

Gauge invariance might be also violated by terms arising 
from momentum dependence of form-factor F. However 
this problem does not arise in the approximation cor- 
responding to equation (|12j) and becomes essential for 
taking into account of g 2 terms. In this case ghost con- 
tributions also do not cancel. The problem of gauge in- 
variance of the next approximations has to be considered 
in future studies. 

We introduce in equation (TT21 an effective cut-off Y, 
which bounds a " low- momentum" region where our non- 
perturbative effects act and consider the equation at in- 
terval [0, Y] under condition 



F(Y) = 0: 



(14) 



and for x > Y we continuously transit to the trivial so- 
lution G — 0. We shall solve equation (fT2|) by iterations. 
That is we expand its terms being proportional to g in 
powers of x and take at first only constant term. Thus 
we have 



F (x) = - -g^afy F Q (y)ydy+-j F (y)dy 
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Fo(y) y 2 dy 
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' r F (y)y^dy)+^^ I F (y)dy. 
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Expression (|15p provides an equation of the type which 
were studied in papers [2H9J , where the way of obtaining 
solutions of equations analogous to (|T5|) are described. 
Indeed, by successive differentiation of Eq.(fT5|) we come 
to Meijer differential equation fill ] 



which solution looks like 

F (x) = ¥„(«) = C x Gll(z |l/2, 1,-1/2, -l) + 
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where 
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is a Meijer function [TI|. In case q = we write only 
indices bi in one line. Constants C\, C2 are defined by 
the following boundary conditions 
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Conditions ([141 IT5|) defines set of parameters 

z = 00 ; C x = ; C 2 = . (19) 

The normalization condition for form-factor F(0) = 1 
here is the following 
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However the first integral in (|2"0"|) diverges due to asymp- 
totic 
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s( z li, 1/2,0, -1/2, -1) l^Z-i 



00 ; 



and we have no consistent solution. In view of this we 
consider the next approximation. We substitute solution 
([TTJl with account of ([20]) into terms of Eq. ([12]) being 
proportional to gauge constant g but the constant ones 
and calculate terms proportional to yfz. Now we have 
bearing in mind the normalization condition 



F{x) = *(z) = l+ 85gA/FVi (lnz + 4 7 + 

96 7T \ 
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where 7 is the Euler constant. We look for solution of 
([2Tjl in the form 

*W = iG?S(z|S il/2>0j _ 1/aj _ 1 )- 

85gVW / 1/2 x 
128vr 15 V z| i,i/2,i/2,-i/2,-iJ + 

CigS2(*I^, 1, -\ -i) + c 2 gS2(*|i, \ -\ -1 

We have also conditions 



1 



z ° , . 87 gVN f Z0 dz 



32 7T 



*(*b) = 0; 



(23) 



and boundary conditions analogous to (118)) . The last 
condition (|23[) means smooth transition from the non- 
trivial solution to trivial one G = 0. Knowing form (|2"2")> 
of a solution we calculate both sides of relation (|2ip in 
two different points in interval < z < zo an d having 
four equations for four parameters solve the set. With 
N = 3 we obtain the following solution, which we use to 
describe QCD case 



g(z ) = 3.8166; z = 0.009553; 
Ci = -5.19055; C 2 = 5.46167. 



(24) 



We would draw attention to the fixed value of parame- 
ter zq. The solution exists only for this value (|2~41 and 
it plays the role of eigenvalue. As a matter of fact from 
the beginning the existence of such eigenvalue is by no 
means evident. This parameter zq defines scale appropri- 
ate to the solution. That is why we take value of running 
coupling g in solution (PHI) just at this point. Note, that 
in what follows we always use the notation F(x) for the 
main form-factor of the approach. 

It is worth to note, that there is also another solution 
of the set of equations, which corresponds to larger value 
of zq ~ 9.6 and smaller value of g(zo) — 0.6 for N = 2. 
We apply this solution for an adequate description of 
non-perturbative contributions to the electro- weak inter- 
action [12Hl5j . 

Let us recall that from three-loop expression for 
a s (/j, 2 ) (JTJ) with number of flavors Nf = 3 we have nor- 
malization of its value at mass of r-lepton © . 

We normalize the running coupling by condition 



a s (x ) 
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1.15515; 



(25) 



where coupling constant g entering in expression (|24j) 
is just corresponding to this normalization point. Now 
from definition of z (fl7| and value zq ([24]) we have 



G = 



A 2 



A G = (2M±7)MeV. 



(26) 



Thus we have obtained the definite value for the coupling 
of the interaction ([5]) under discussion. Typical energy 



scale around 250 MeV is natural for strong interaction. 
It is also worth mentioning the value of the momentum 
which corresponds to boundary of non-perturbative re- 
gion zq. From Eqs. (|24[ I26p we have for this momentum 



(27) 



po = (630 ± 18) MeV . 



Non-perturbative boundary (|27p seems also natural from 
phcnomcnological point of view. 

We have to bear in mind, of course, that all these re- 
sults are obtained under chosen approximation. For ex- 
ample, change of form of dependence on three variables 
in expression (|10p leads to some change in constant term 
in inhomogeneous part of equation (|2"Tj) . The coefficient 
afore the logarithm in its second term does not depend 
on the form, but the constant one can be changed. It is 
important to understand how small changes in this term 
influence results. In view of this we consider additional 
term e in the inhomogeneous part of (|2ip . Thus we have 
the following modified expression 
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Let us take example e = 0.13. In this case instead 
of (HJU) we have 

g(z ) = 3.11587 ; z = 0.0153348; 

Cx = - 4.47289 ; C 2 = 3.62922 ; (29) 

that in the same way as for case e = leads to the 
following parameters 

a s (xo) = = 0.7726; 

4-7T 



G = -i- ; A G = (273.5 ± 7.0) MeV . (30) 



Another example e — 0.15. In this case we have 

g(z ) = 3.03685; z = 0.0163105; a s (x ) = 0.7339; 
Ci = - 4.37005 ; C 2 = 3.43372 ; 

G = - T ; A G = (276.4 ±7.0) MeV. (31) 



III. RUNNING COUPLING 

In previous sections N.N. Bogoliubov compensation 
principle Q was applied to studies of a spontaneous 
generation of effective non-local interaction © in QCD. 

It is of the utmost interest to study an influence of 
interaction ^ on the behavior of strong running coupling 
a s (k 2 ) in the region below zq i.e. k < Pn (j2"7| . 
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FIG. 2: Diagrams, describing the contribution of non- 
perturbative vertex denoted by the black spot, to the 

running coupling a B (k 2 ). Simple lines correspond to gluons 
and thick lines correspond to quarks. 



For the purpose we rely on considerations connected 
with the renormalization grou p approach [l6[ (for ap- 
plication to QCD see, e.g [17| )■ We have the one loop 
perturbative expression for QCD ^-function. 



(4 7T) 
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(32) 



We shall take additional contributions for small momen- 
tum k 2 — > of our new interactions according to dia- 
grams shown in Fig. 2 that gives instead of 



(4tt) 
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2N f \ 405 v / 3 5 (z0) 



2tt 



*(0) 



; (33) 



where $(0) is the result of calculation of diagrams Fig. 2 
(see below) . Here we see a decisive difference in behavior 
of perturbative f3 (|3"2"1) . which acts at large momenta 
k > po and non-perturbative one for small k ~ (|33|) . 
According to calculation of $(0) with account of (pM|) the 
sign of (3 changes between these regions. So a s (k 2 ) for 
k 2 — > is also positive as well as for large k. To consider 
a behavior in between we return to definition of the /3- 
function (171 
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where function $(i) is defined by calculation of diagrams 
Fig. 2 and d — > 4 is the space-time dimension. 



*(*) = 



u-U/A 



t u-t/2 
1 4(u-3t/4) 2 



F(u)du + 



3t/4 t(u-t/2) 

zm 4(u-3i/4) 2 



F(u)du; t < zOl; 



(35) 



= 0; i > 



3i/4 t(u-t/2) 

—5— ; zoi — v z o : 



F(u)du; zqi < t < 



4z i 



%/3Gq 2 



%/3Gfc 2 



32 7T ' 32 7T 

This leads to modification of relation (|33p 



P(g,t) = - 
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(47T) 



fn 27V; 405V3g(zO) 



2tt 



$(t) 



; (36) 



Thus in approximation using the two-loop expression 
corresponding to diagrams of Fig. 2 we have for Nf = 3 



a s (k 2 ) = 



«(*o) 



1 + 



6tt ^ 



21 

2 



2tt 



Ht) In 1 



(37) 



where t is defined in (|3"5j) . With G defined by (|2^|) . 3(^0) 
defined by ^ and fc 2 = Q 2 we have the behavior of 
a s (Q). With fixed parameter e in (12^1) we calculate the 
behavior of running coupling. Let us begin with initial 
case e = 0. We have value of a s (Q) at the beginning 
point of non-perturbative contribution, corresponding to 
zqi = |zoij corresponding to momentum Q — 726 MeV. 



a s {z 01 ) = 0.936. 



(38) 



The boundary of non-perturbative region Qo = 726 MeV 
seem quite reasonable. 

Now the behavior of a s (Q) is drawn in Fig. 3. We 
would like to draw attention to the result, presented at 
Fig. 3 , which consists in absence of Landau pole in ex- 
pression (|37l) . Remind, that in perturbative calculation 
up to four loops the singularity at Landau pole point 
is always present. Only by taking into account of the 
non-perturbative effects we achieve elimination of this 
very unpleasant feature, which was seriously considered 
as a sign of the inconsistency of the quantum field the- 
ory [H,|2|. 

There is also a feature of expression (13"T|) . which de- 
serves being mentioned. The limit of a s (Q) for Q — > 
is zero. Such possibility is also discussed on phenomeno- 
logical grounds. In particular, there are indications for 
decreasing of a s (Q 2 ) for Q 2 —> in studies of low-mass 
resonances [Isj . Let us note, that a number of lattice 
calculations of the running coupling also give a similar 
behavior [lj-[23|. 

Let us also consider a s behavior for other values of 
parameter e. The behavior for e = 0.13 is presented in 
Fig. 4. The pole here is also absent, but values of a s in 




FIG. 3: Dependence of the running coupling a s (Q), Q in 
MeV, with e = 0. The continuous line corresponds to a s 
with non-perturbative contribution (|37[) . the discontinuous 
one with a pole corresponds to the usual perturbative one- 
loop expression. 
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FIG. 5: Diagrams for calculation of the gluon condensate. 
Lines - gluons, black circle - triple vertex (J7]), open circle 
- four gluon vertex with corresponding form-factor and 
skew cross - vertex (|4ip . Momenta directed to the right are 
p-q/2, q, -p-q/2 for bug-like diagrams and p-q/2, p+q/2 for 
oo-like diagrams. 



of operator (f^U|) in brackets. We denote it by skew cross 
in Fig. 5 



(41) 



With distribution of integration momenta denoted in 
Fig. 5 form-factor in both types of diagrams according 
to (|10p has the same argument: 



F( P 2 



(42) 



FIG. 4: Dependence of the running coupling a 3 (Q) for 
e = 0.13. The continuous line corresponds to a s with non- 
perturbative contribution (|37p . the discontinuous one with a 
pole corresponds to the usual perturbative one-loop expres- 



non-perturbative region are smaller than for case e = 0. 
The average a s in the non-perturbative region for e = 
0.13 

a s = i- [ Q ° a s {Q)dQ = 0.87; 

VO Jo 



'1287TZ01 



(39) 



Y 3V3G ' 
For e = 0.15 a s = 0.84. 

IV. THE GLUON CONDENSATE 

One of important non-perturbative parameters is the 
gluon condensate, that is the following vacuum average 



V 9 = < -9— F a F a > 



(40) 



Let us estimate this parameter in our approach. We ap- 
ply our method to the first non-perturbative contribu- 
tions, presented at Fig. 5, which is proportional to gG. It 
is important to introduce Feynman rule for contribution 



It comes out, that the second and the third terms in 
the second row of Fig. [5] are twice each of the previous 
terms. Thus the sum is equal to the result for the first 
diagram multiplied by 10. 

We have after the Wick rotation 



Vn 



10 x 24 g 3 G 

(2 7T) 8 7T 2 



(43) 



F( P 2 



79 



12 (p 2 q 2 



pq 



4" ' q 2 {p- q/2) 2 (p + q/2) 2 
Using the following integral by angle 
sin 2 {6)d6 



dpdq . 



o (p 2 + q -) 2 



2(z+§) 



- (pq) 2 



(44) 



>J 



, +e( 

4 7 x V 4 



y 



x =p ; y=q; 

we obtain the following expression for quantity (|43[) 
5«7 3 2 n 



V 2 



It = 12 - 



G 2 7T 3 y/S Jo 

l (t-y) 2 



o t-y/2 

t (t - y/ 2 )y 

GV3( 3y 

t = : X + 

2 5 7T I 4 



F{t) I t dt 
dy - 
dy + 



it/3 



(45) 
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We have already expressions (|22[ 124)) for form-factor 
F(z), z — t 1 . So calculation here is direct and we obtain, 
using values for g (|24|) and the central value in definition 
of G (H| 



Vn 



5g 3 2 10 
tt 3 V3G 2 



0.00955 GeV 4 



u(2 -61n^ J F{z)^fz~dz 



(46) 



Provided we take nonzero value for e in expression ([28 
results for gluon condensate read 



V 2 = 0.0120 GeV 4 (e = 0.13); 
V 2 = 0.0128 GeV 4 (e = 0.15) 



(47) 



So in this approximation we have the non-zero non- 
perturbative parameter Vi. Its value agrees within ac- 
curacy of determination of this parameter with phe- 
nomenological values V 2 ^ 0.012 GeV 4 0, V 2 ^ 
0.010 GeV 4 [H]. Values flU g7]) show variation in the 
range of uncertainty of its phenomenological definition. 
Thus we can state, that our non-perturbative approach 
allows to calculate safely this important parameter. 
Let us also estimate vacuum average V3 

V 3 = < g 3 f abc F£ v F h vp F^ > . (48) 

Quite analogous calculations give e.g. with e = 0.13 



o 3 2 17 / 4\ f z ° 

V 3 = ^3- (2 -6 In -J J zF(z)dz 

0.00744 GeV 6 . 



THE GLUEBALL 



(49) 



The existence of anomalous interaction ([6J makes pos- 
sible to consider gluonic states. We shall consider scalar 
glueball Xq state to get indications if value of the non- 
perturbative constant (1261) may be used for adequate de- 
scription of the non-perturbative effects of the strong 
interaction. For the purpose we use Bethe-Salpeter 
equation with the kernel corresponding to one-gluon ex- 
change with our (point-like) anomalous three-gluon in- 
teraction ©. We take for vertex of Xq interaction with 
two gluons in the following form 

^fF^F^XoV^x); x=p 2 - (50) 

where <f g b(x) is a Bethe-Salpeter wave function. We have 
for the first approximation (zero momentum of X ) 



3G 2 / 1 



16tH \2x 2 J 



-[ y 2 ^ gb (y)dy-3 [ y^ gb (y)dy- (51) 
x Jo Jo 

X [ Y ^y)dy + ^[ Y ^dy, 
Jx z Jx v 



where we take again the upper limit Y of integration as 
in (|T2l due to form-factor of interaction © F(x) = for 
x > Y. Again by successive differentiations we obtain 
from Eq. (f5"T]) the following differential equation 



dz' \ dz' 2 A dz' 2 



\z'^-l\y gb {z') = z'y gb {z') + 



(52) 



f z ° 9G 2 r 2 9G 2 u 2 

G = 4/ * gb (t')dt'; z' = f = JUV 

Jo 



1287T 2 



1287T 2 



Comparing variable z' in Eq. (|52p with the initial variable 
z in Ea dTTj) we see relation z' — 24 z. This means also, 
that zo = 2Azq, zq from solution (|2"4"1) . In new variables 
Eq. (l51[) . in which we also have taken into account terms, 
proportional to gauge coupling g and mass of the bound 
state squared to 2 , looks like 



= 1 - / * 9 b(t')t'dt' 

^ z JO 



3 J z , V 
1 = 4^ \ gb {t')dt' + (k 



3gV2 
2tt 



(53) 
^Idt'. 



Here k is connected with the bound state mass to in the 
following way 



3Gto 2 
8V2tt 



(54) 



According to expression (|52|) we look for the solution of 
Eq. (l5"Tj) in the following form 

= £^5 0,1/2,-1/2,-1) + (55) 
CxGf^Al, l/2,-l/2,-l) + 
C 2 GH(-z'\l, 1/2, -1/2, -l) . 



By substituting expression (I5"5"j) into set of equations ([53")) 
and using the values of g and zq (|24p we obtain unique 
solution for parameters 

Gi = 1.07899; G 2 = -1.38099; k = -2.6415. (56) 



Now from values ([261 156]) . using relation (|54l) . we have 
the lightest scalar glueball mass 



to = 1479 ±40MeF. 



(57) 



This value is quite natural, the more so, that the most 
serious candidate for being the lightest scalar glueball 
is the state Jo (1500) (see recent review (26J) with mass 
1507 ± 5MeV, that evidently agrees our number ([5"T]) . 
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Now we have to obtain the coupling constant of the 
scalar gluon entering in the expression of the effective in- 
teraction (|50|) . For the purpose we use the normalization 
condition for Bethc-Salpeter wave function ^(t). 



1 



dt' . 



(58) 



Substituting into Eg. ((58)) solution ([55 
ing the integral, we obtain 



and calculat- 



r 2 

I = 



V2I ' 



1.825 G ; 
dt' = 1.21732. (59) 



From result (I5"9l we have the following value of the glue- 
ball coupling 



G gb - 



1 



190.337 MeV 



5.254 
GeV 



(60) 



VI. CONCLUSION 



An existence of a non-trivial solution of a compen- 
sation equation is extremely restrictive. In the most 



cases such solutions do not exist at all. When we start 
from a renormalizable theory we have arbitrary value for 
its coupling constant. Provided there exists stable non- 
trivial solution of a compensation equation the coupling 
is fixed as well as the parameters of this non-trivial so- 
lution. Note, that application of the same approach to 
the electro-weak theory [L2T - [l5l . |27| also leads to strong 
restrictions on parameters of the theory including the 
coupling constant. 

We also may state, that in the case, discussed in the 
present paper, just the non-trivial solution is the stable 
one, because the theory with the Landau pole is unstable. 

We consider the results for the gluon conden- 
sate (146I47|) and the glueball mass (|57p as a confirma- 
tion of efficiency of our approach in application to non- 
perturbative contributions to QCD. 

Thus we consider the present results for low-momenta 
a s to be encouraging and promising for further appli- 
cations of the Bogoliubov compensation approach to 
principal problems of elementary particles physics. 
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